2011 Fermi Symposium, Roma., May. 9-12 



A variability study of the AGILE first catalog of 7-ray sources on 2.3 
years of AGILE pointed observations 

F. Verrecchia^'^, C. Pittori^-^ A. Bulgarelli*, A.W. Chen^, M. Tavani*5, F. Lucarelli^-^, P. Giommi\ on 
behalf of the AGILE Collaboration 

1 ASI Science Data Center, ESRIN, 1-00044 Frascati (RM), Italy 

2 Consorzio Interuniversitario Fisica Spaziale, villa Gualino - v.le S. Severe 63, 1-10133 Torino, Italy; 

3 INAF-OAR, Astronomical Observatory of Rome, Monte Porzio Catone, Italy; 

4 INAF-IASF Bologna, via Gobetti 101, 1-40129 Bologna, Italy; 

5 INAF-IASF Milano, via E. Bassini 15, 1-20133 Milano, Italy; and 

6 INAF-IASF Roma, via del Fosso del Cavaliere 100, 1-00133 Roma, Italy; 



AGILE pointed observations performed from July 9, 2007 to October 30, 2009 cover a very large time interval, 
with a 7-ray dataset useful to perform studies of medium to high brightness galactic sources in the 30 MeV 
- 50 GeV energy range. We present a study of the lAGL galactic sources in the E > 100 MeV band, over the 
complete Agile pointed Observation Blocks (OBs) archive. The first AGILE Gamma-Ray Imaging Detector 
(GRID) catalog included a sample of 47 sources (lAGL; |19|I ). detected with a conservative analysis over the 
first year of operations dataset. In the analysis here reported we used data obtained with an improved full Field 
of View (FOV) event filter, on a much larger (about 27.5 months) observation dataset, analyzing the merging 
of all data and each OB separately. The data processing resulted in an improved source list as compared to the 
lAGL one, particularly in complex regions of the galactic plane. We present here some results on the revised 
lAGL galactic sources and on the variability of some of them. 



1. Introduction 

AGILE (Astrorivelatore Gamma ad Immagini LEg- 
gero) ([22,]) is a mission of the Italian Space Agency 
dedicated to 7-ray astrophysics in the 30 MeV - 50 
GeV energy band, and simultaneously to hard X-ray 
in the 18 - 60 KeV band. AGILE, in orbit since April 
23 2007, has been the first instrument of a new gen- 
eration of high-energy space missions based on the 
solid-state silicon technology, permitting to advance 
our knowledge on many source classes from active 
galactic nuclei, to pulsars (PSRs), unidentified 7-ray 
sources, Galactic compact objects and supernova rem- 
nants. On June 11, 2008 the Fermi Gamma- Ray Space 
Telescope ([13,01) 

was launched, and is currently op- 
erating together with AGILE. 

The AGILE Payload detector consists of the silicon 
tracker (ST; [1, jl^), the X-ray detector SuperAG- 
ILE (0), the CsI(Tl) Mini-CalorimeterfMCAL; Q) 
and an anti-coincidence system (ACS; The com- 

bination of ST, MCAL, and ACS forms the Gamma- 
Ray Imaging Detector (GRID). GRID is sensitive to 
photon energies in 30 MeV - 50 GeV energy band, 
and thanks to silicon technology has a wide FOV (2.5 
sr in pointing) and accurate timing (a few /xs) and po- 
sitional information (15' location accuracy for > lOcr 
detection), with an angular resolution of 3.5° at 100 
MeV and 1 ° above 1 GeV. 

AGILE spacecraft operated in fixed-pointing 
mode since October 2009 (completing 101 pointings 
or Observation Blocks, OBs), when the attitude 
control system had to be reconfigured into "spinning 
operation mode" . Currently the instrument pointing 
direction scans the sky with an angular velocity of 



about 0.8°/s, accessing about 80% of it each day. The 
AGILE satellite raw data, down-linked about every 
100 minutes, are transmitted from Mission Control 
Center at Telespazio, Fucino, to the AGILE Data 
Center (ADC), part of the ASI Science Data Center 
(ASDC) located in Frascati (Italy). The ADC has 
the duties of data reduction, scientific processing and 
archiving and finally to distribute standard Level-2 
data to Guest Observers (GOs) or, when data become 
public, to all the scientific community (see ADC web 
page http://agile.asdc.asi.it) . 

The significance-limited (4 sigma) sample of 47 
lAGL sources were detected in the E > 100 MeV band 
with a conservative analysis of the inhomogeneous 
first-year sky coverage dataset. We present here the 
first results of a variability study of a sample of 62 
sources (Fig. 1) analyzing separately each OB in the 
2.3 years AGILE pointing mode dataset. The sam- 
ple was obtained with a revision of the lAGL sample 
on the maps obtained from the whole dataset ( "deep" 
maps from now on). 



2. OB source detection procedure 

The standard analysis OB pipeline at the ADC 
was executed at the end of an OB to remove data 
corresponding to slews and occasional losses of fine- 
pointing attitude, and to build the official OB data 
archive. Moreover scientific maps in E> 100 MeV en- 
ergy band with a size of 60° x 60°, binned at 0.25°, 
were created selecting confirmed events and exclud- 
ing albedo contaminated time intervals. The "point- 
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orange: 

red: 8<SQRTTS. 
green: 5<SQRTTS<8 
blue: SQRTTS<5 



Figure 1: The AitofT plot of the 62 distinct source positions detected in a revision of the original 47 lAGL ones, on all 
the pointed observations data (colors and symbol sizes are proportional to significance, blue the lowest and orange the 
highest). On the right the AGILE-GRID 2.3 years all-sky exposure map in Aitoff projection from all the pointed 
observations. 



ing mode" OB archive is composed of 101 OB cov- 
ering the wide timespan with non uniform exposures 
(ranging from Id to 45dd). All the OB were recently 
reprocessed with the last software release. The proce- 
dure developed for this analysis on the whole archive 
in E> 100 MeV band, is based on source detection 
at fixed preselected positions using the AGILE Maxi- 
mum Likelihood (ML). It consists of two main steps: 
I) preliminary revision of the lAGL source list, based 
on updated deep maps and on more recent published 
AGILE results; II) execution of a ML multi-source 
task on each OB data in a specific pipeline, keep- 
ing reference source positions fixed to those previously 
determined. In the analysis no automatic source de- 
tection process is implemented. A check leaving the 
position free is implemented to verify the goodness of 
candidate transient sources. Both steps were repeated 
for 5 main iterations to improve the multi-source de- 
tection in complex regions (see section 12. 2p . 



2.1. Source detection method 



The detection method used is the ML. The signifi- 
cance (number of sigmas) of a source detection is given 
by the square root of the test statistic TS, defined as 
-2 times the log of the likelihood ratio, and expected 
to behave as ([3, Chen et al. in preparation). 
The likelihood ratio test is built considering for the 
background only hypothesis the AGILE diffuse 7-ray 
model ([m, Giuliani et al. in preparation). We used 
the task "AG_multi2" included in the AGILE software 
for the ML "multi-source" analysis. 



2.2. The First AGILE catalog revision on 
the 2.3 yrs dataset 

The first AGILE catalog ([l^) was built from data 
of the first year (July 2007- June 2008). The data 
analysis, based on a conservative event filter ([16]), 
and the non uniform sensitivity due to inhomoge- 
neous first-year AGILE sky coverage, limited the re- 
sults in complex galactic regions. A preliminary revi- 
sion of lAGL sources in regions such as the Carina, 
Cygnus, Crux and Galactic Center fields, was realized 
on updated merged maps from data up to October 
2009 (see Fig. [1]), using data from a new event fil- 
ter (FM3.119; Bulgarelli et al in preparation). In this 
revision new candidate source were considered, and 
some with 3 < ^/TS < 4 were also included in the ref- 
erence list to check their variability among all the OBs 
(see Fig. ^. 

Moreover recent AGILE results obtained in these 
regions were taken into account ([23|, [13], [lH, [l3l |. 
[1]). We considered also preliminary results from the 
ADC Quick Look processing on week timescale and 
the procedure results in the first 4 iterations. We ob- 
tained so an updated reference source list. 

A new source detection procedure on deep maps of 
all AGILE data in pointing mode has been realized for 
the next AGILE catalog and results will be reported 
in Bulgarelh et al. 2011. 

3. Detection selection and variability 
analysis 

The results from the last execution (the fifth iter- 
ation) of the ML multi-source task on sources in the 
reference list, were filtered taking into account off-axis 
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AGILE GRID 2007 Jul 9 Exp'jsure; 71J 



A5ILE ^RID 2007 Jul 9 Exposure: 715S49O0 s 





Figure 2: On the left a simple case of position refinement, the high latitude lAGL J1846+6714. The new smaller error 
circle is indicated inside the lAGL one. Right, a difficult position refinement case, the Carina region . New positions 
are indicated with small squares, with sizes not proportional to errors, while the lAGL error circles are reported. 



Crab, OB Processing 




5.44x10* 5.46x10'' 5.48x10' 5.5x10* 

TIME (MJD) 

Figure 3: The Crab E>100 light curve (in 10"* 
ph cm~'^ s~^; upper panel) and the plot of y/TS and the 
off-axis angle (in degrees) vs time (in MJD). The lAGL 
flux value and the 1 sigma error levels are shown as 
magenta and black dashed lines. The 2007 flare episode 
is the most significant detection, at a fiux of 619 ± 71 
ph cm~^ s~^. 



angle and significance. Our goal was to investigate 
the variability within the OB dataset for established 
sources on deep maps, and possibly search for flare 
episodes for some low significance source. As final 
selection only those sources having at least one de- 
tection at significance higher than Su were accepted, 
for all sources having high significance (>=4) on deep 
maps, or otherwise with at least one detection above 
4cr. Then all detections with y/{TS) > 2 were con- 
sidered in the variability analysis. Source 7-ray flux 



variability among all OBs was tested according to the 
method developed by [l^ and recently reported in an 
analysis of lAGL J2022-h4032 data ([1|). The vari- 
ability index V=-log(Q) is evaluated after the calcu- 
lation of the weighted mean flux and its error and so 
also the relative x^. Q is the probability of having a 
value of > Xobserved ^'^^ ^ sourcc with constant flux. 
Sources with V<0.5 are usually classified as "non vari- 
able", with 0.5<V<1 as "uncertain" and with V>1 
as "variable" . Another index was calculated for com- 
parison (used in Fermi catalogs; Q), as a simple 
■ Both indices were computed also adding a system- 
atic component of 10% to flux errors {Vsys, x'^sys), to 
take into account the not well known particle back- 
ground systematic, expected to vary among different 
integrations. 



4. Conclusions and future developments 

Applying our selection criteria, we obtained a sam- 
ple of 1267 detection of 62 distinct sources (Fig. [TJ 
previous processing preliminary results were presented 
in [2^). For all sources a revision of single detec- 
tions is on-going as a verification and improvement of 
the multi-source ML analysis and taking into account 
more recent calibrations. The non uniform exposures 
among the OBs, from 1 to 45 days, together with the 
pointing strategy, put strong constrains to the capac- 
ity of "resolving" complex regions in a large number 
of OBs and also to the variability analysis. 

The light curve of the most bright sources with sta- 
ble flux, such as Geminga and Vela PSR, have been 
preliminarily checked. A particular case is the one of 
the Crab, which has been recently reported to have 
a "non stable" 7-ray emission ([5^, |26j, Q, [30l|). 
In Fig. [3] is shown the Crab light curve showing the 
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October 2007 flare episode, discussed in [26|, at an 
higher flux compared to the 2010 flare. We obtained 
so a \ sys of 0.4 for Geminga and 3.36 for Crab. 

As an example of long-term monitoring light curves 
in the E>100 MeV band obtained in this analysis 
at the OB timescale, we show below (Fig. 5 and 6) 
the results for the two known galactic sources lAGL 
J1836+5923, associated to the LAT PSR J1836+5925, 
and lAGL J0242+611 associated with the HMXB LS 
1+61303. 



1AGLJ1836+5923, OB Processing 



The use of Ysvs and 



gave compatible results 
for most of the sources. The source class which in- 
cluded more variable sources were the bright Blazars 
according to both parameters, the most variable one 
being 3C 454.3. The complete refined results on the 
OB timescale will be presented in Verrecchia et al. 
2011. 
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Figure 4: the lAGL J1836+5923 and lAGL J0242+611 
light curves in E>100 MeV band, together with the y/TS 
and the ofF-axis angle plots. The magenta line is the 
lAGL flux and the two black dashed lines are the levels 
for 1 o" errors 
respectively. 
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